Assessing the Polarization of a Quantum Field from Stokes Fluctuations 
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We propose an operational degree of polarization in terms of the variance of the projected Stokes vector 
minimized over all the directions of the Poincare sphere. We examine the properties of this degree and show 
that some problems associated with the standard definition are avoided. The new degree of polarization is 
experimentally determined using two examples: a bright squeezed state and a quadrature squeezed vacuum. 
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Introduction. — Polarization is a fundamental property of 
light that has received a lot of attention over the years IH]. 
Nowadays, the topic is witnessing a revival in interest because 
of the fast developments, both on applications and on more 
fundamental aspects. As polarization is a robust character- 
istic, relatively simple to manipulate without inducing more 
than marginal losses, it is not surprising that many experi- 
ments at the forefront of quantum optics involve this observ- 
able n. 

In classical optics, polarization can be elegantly visualized 
by using the Poincare sphere and is determined by the Stokes 
parameters. These are directly measurable quantities that can 
be straightforwardly extended to the quantum domain, where 
they become operators |3fl . 

The classical degree of polarization is simply the length of 
the Stokes vector. While this provides a very intuitive pic- 
ture, for many complex fields it has also serious drawbacks. 
Indeed, this classical quantity does not distinguish between 
states having remarkably different polarization properties 10] • 
In particular, it can be zero for light that cannot be regarded 
as unpolarized, giving rise to the so-called "hidden polariza- 
tion" n. 

These flaws have prompted some novel generalizations of 
the degree of polarization A notion that has been 

gaining support in the quantuni optics community is to ap- 
ply a properly chosen distance lll2n (entropy can be regarded 
as a special case lO ). This has the potential advantage of 
circumventing most of the aforementioned difficulties, while 
making close contact with other measures introduced to quan- 
tify quantum resources HJl . 

There is, however, a problem with this approach: these dis- 
tances can be computed (not measured) only after a complete 
knowledge of the state, which in practice implies a full quan- 
tum tomography. In other words, while offering very good 
properties, distance measures do not have a clear operational 
meaning. 

We adhere to the view that the Stokes variables constitute 
a natural tool in appraising polarization properties, so they 



should be the basic building blocks for any practical degree 
of polarization. One can expect that the problems arising with 
the classical degree are due to its definition in terms exclu- 
sively of first-order moments of the Stokes variables. This 
may be sufficient for most classical problems, but for quan- 
tum fields higher-order correlations are crucial. 

Our goal in this Letter is to provide such a characterization. 
We learn from coherence theory that a full description of in- 
terference phenomena may involve a hierarchy of degrees. In 
this vein, we go beyond the first-order description and look for 
a second-order degree as the minimum Stokes variance over 
all directions of the Poincare sphere. This simple proposal 
will prove very satisfactory when facing the complications 
known in this field. We also present a couple of experimental 
examples confirming the feasibility of our scheme. 

Polarization structure of quantum fields. — We begin by 
briefly recalling some background material. We assume a 
two-mode quantum field that is described by two complex am- 
plitudes, glh and ay, where the subscripts H and V indicate 
horizontally and vertically polarization modes, respectively. 
The commutation relations of these operators are standard: 
[aj,a|.] = 8jk, with j, /c G {H^V}. The analysis is greatly 
simplified if we use the Stokes operators 
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together with the total photon number So = N 
aydy- The average values of these operators are precisely the 
classical Stokes parameters. One immediately finds that the 
components of the Stokes vector S = {Sxj Sy^ SzY (where t 
denotes the transpose) satisfy the commutation relations dis- 
tinctive of the su(2) algebra: [Sx^ Sy] = i2Sz and cyclic per- 
mutations. This noncommutability precludes their simultane- 
ous precise measurement, which is expressed by the uncer- 
tainty relation 

{ASf = {ASxf + {ASyf + {ASzf > 2{So) , (2) 
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with (ASk)'^ (k = x^y^ z) being the corresponding variances. 
In addition, [^o, S] = 0, so we can treat each subspace with 
a fixed number of photons N separately. This can be em- 
phasized if instead of the Fock basis for both polarization 
modes, \n)H\'m)v (n^m = 0, . . . , oo), we employ the rela- 
beling \N,k) = \k)H\N-k)v (/c = 0,1,...,7V). In this 
way, for each fixed A^, these states span an SU(2) invariant 
subspace of dimension N 

The standard definition of the degree of polarization for a 
quantum state g is 



(So) 



(3) 



where the subscript 1 stresses here that it involves first-order 
moments of the Stokes variables. We note that for any single- 
mode state of the form |^)i:/|0)y, we get Pi = 1, which 
seems unphysical for a variety of reasons. In particular, when 
1^)// |0)//, we have Pi = 1 for field states arbitrarily close 
to the quantum two-mode vacuum. Moreover, unpolarized 
states according to © are determined by (S) = 0. Nonethe- 
less, there are states fulfilling this latter condition (as, e.g., 
1^)// ) that cannot be regarded as unpolarized, as revealed 
by a number of features. These unwanted consequences call 
for alternative measures. 

Second-order quantum degree of polarization. — As from 
the previous discussion, it seems clear that higher-order 
moments must be taken into account, as advocated by 
Klyshko fisl]. For the time being, we concentrate on the sec- 
ond order: our task is thus to link the resulting fluctuations 
with our notion of a polarization degree. To this end, we ob- 
serve that a sensible modification of ([3]) is easily obtained by 
replacing (^o) with [(^o(^o + = (S^) V2 in the de- 

nominator. The resulting degree llldl 



{ASl 
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contains the desired second-order information (hence the sub- 
script 2) and fixes some of the above-mentioned problems. 
For example, P2 < 1 for every state |^)i:/|0)y, and P2 ^ 
when \'^)h \^)h- However, other bugs still persist. The 
reason is that Q does not properly represent the behavior of 
the fluctuations in phase space. To catch these aspects we pro- 
pose to use 




inf 
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where Sn = S - n, with n being a unit vector in an arbitrary 
direction of spherical angles (6>, (j)). The factor 1/3 has been 
introduced for normalization. 

To further appreciate this idea, we define the real sym- 
metric 3x3 covariance matrix for the Stokes variables as 
= UiSk.SA) - {Sk){Si), where {,} is the anticom- 
mutator llTII . In consequence, (A5'n)^ = n^Fn and, since 



r is positive definite, the minimum of (A5'n)^ exists and it 
is unique. If we incorporate the constraint n^n = 1 as a La- 
grange multiplier 7, this minimum is given by Fn = 711: the 
admissible values of 7 are thus the eigenvalues of F and the 
directions minimizing (A/Sn)^ are the corresponding eigen- 
vectors, which, following the standard nomenclature in statis- 
tics, are known as principal components of F. 

The covariance matrix F can be made diagonal by an or- 
thogonal matrix R. In this rotated reference frame we have 

that «S = RS satisfies «S =8^, so that 

{ASif + + (ASsf < {S') = (S^) , (6) 

where the subscripts 1, 2, and 3 indicate the directions of the 
orthogonal eigenvectors of F. The contour surface of these 
variances defines an ellipsoid that provides an accurate repre- 
sentation of the noise distribution of the state. 

Properties and examples. — Let us explore some proper- 
ties of the degree P2 . Unpolarized states according to P2 are 
those whose fluctuations are isotropic and saturate the bound 
in Eq. Q. This means that the ellipsoid reduces to a sphere 
of a radius (KS^))^/^ As we shall see, this allows to dis- 
tinguish hidden polarization not revealed by Pi. We note, in 
passing, that the unpolarized states introduced in Ref. as 
those invariant under SU(2) transformations are also unpolar- 
ized for P2. However, the converse is not true in general. 

It follows directly from the definition of P2 that any SU(2) 
polarization transformation U leaves the degree of polariza- 
tion invariant: ¥2(0) =F2{U gU"^). 

It is clear that the moments of any energy-preserving ob- 
servable (such as S) do not depend on the coherences between 
different subspaces. This means that the only accessible infor- 
mation from any state g is just its polarization sector, which is 
defined by the block-diagonal form ^poi = X]^=o ^ a'^ ^ ^ a'^' 
where Iat is the projector onto the TV-photon subspace. 
Therefore, any g and its associated block-diagonal form ^poi 
have the same degree of polarization P2. This is consistent 
with the fact that polarization and intensity are, in principle, 
independent concepts: in classical optics the form of the el- 
lipse described by the electric field (polarization) does not de- 
pend on its size (intensity). All this confirms that our proposal 
fulfills all the requirements for a bona fide second-order de- 
gree of polarization. 

We further develop these ideas by presenting a few relevant 
examples. First, for any two-mode number state \n)H\'m)v, 
^2{\n)H\'m)v) = 1- In particular, this means that P2 identi- 
fies the hidden polarization of, e.g., the state |n)if |n)y. 

For two-mode quadrature coherent states \a)H\P)v, with 
an average number of photons TV = |ap + simple cal- 
culations give F2{\a)H\f^)v) = [TV/(TV_+ 3)]^/^ so when 
TV 00, P2 tends to unity, and when TV 0, P2 tends to 
zero, showing a good classical limit. Interestingly, the covari- 
ance for these states is isotropic and the corresponding ellip- 
soid reduces to a sphere of radius TV^/^. 
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FIG. 1 : (Color online) (Top) Setup for efficient generation of a po- 
larization squeezed state and the corresponding Stokes measurement 
apparatus. (Bottom) Measured variances for that state with the indi- 
cated scale (in dB noise power relative to the shot noise, marked by a 
white line). The minimum measured variance is — 5.0 db 0.3 dB. The 
white point is the tip of (S). We include also a zoom around the Sx 
axis, near the minimum variance. 



FIG. 2: (Color online) (Top) Setup for the generation of a squeezed 
vacuum in the horizontal mode and vacuum in the vertical one. We 
measure -3.8 dB (8.6 dB) of quadrature squeezing (antisqueezing) in 
the horizontal mode. (Bottom) Measured variances {/\SnY for the 
state as a color map on the Poincare sphere with the indicated (linear) 
scale. 



For single-mode states \^)h\^)v^^^ find 



P2(|^)/f|0)x 



3min[(A7V)2,A/'] 
(A7V)2 +7V(7V + 2) 



(7) 



where TV is the average number of photons. The problems 
arising with these states when using Pi are thus avoided. 

We finally consider SU(2) coherent states |^,(/ 
D{e,(t))\N,k = 0), where £>((9,0) = exp[6>/2(6'+e- 
S-e'^^)], with S± = Sx ^ iSy, is the standard displacement 
operator on the sphere. These are the only ones that satisfy 
relation dJ]) as an equality, so P2(|^, = 1 and they are 
completely polarized for our approach. Incidentally, we have 
also Pi (|6>, 0)) = 1 and one could expect that they would ful- 
fill similar properties for arbitrary orders. 

Experiment. — We demonstrate our proposal with two dif- 
ferent quantum states of light: a very bright polarization 
squeezed state and a quadrature squeezed vacuum produced in 
an optical fiber and in an optical parametric oscillator (OPO), 
respectively. 

To generate the bright squeezed light, we employ ultra- 
short laser pulses in the soliton regime of an optical fiber to 
achieve a large effective nonlinear Kerr response and avoid 
dispersive pulse broadening (see Fig. 1) Our experiment 
uses a Cr^+:YAG laser emitting near Fourier-limited 140 fs 
FWHM pulses at 1497 nm with a repetition rate of 163 MHz. 
We utilize the two polarization axes of a 13.2 m birefringent 
fiber to simultaneously generate two independent quadrature 
squeezed states in the H and V modes, with a relative phase 



of 7r/2. The average output power from the fiber was 13 mW, 
which, with the bandwidth definition of our quantum state, 
corresponds to an average number of photons of 10^^ per 1/is. 

The Stokes measurement is also shown in Fig. 1 and con- 
sists of a half- wave plate (A/2, 9) followed by a quarter- wave 
plate (A/4, (j)) and a polarizing beam splitter (PBS). The trans- 
formation performed by the wave plates can be represented by 
l)(6>, 0), while the PBS projects on the basis \N^k). The out- 
puts of the PBS are measured using high efficiency photodi- 
odes (98%), the photocurrent difference is produced, and the 
resulting fluctuations are evaluated at a sideband of 17.5 MHz 
(and a bandwidth of 1 MHz). In this way, the setup enables 
the measurement of S^i lEoll . 

For each pair of angles (6>, (j)) the noise statistics were ac- 
quired and the optical intensities incident at both detectors 
were recorded. From this, the Stokes variances (A5'n)^ were 
obtained and the results are plotted in Fig. 1 as a color map on 
the sphere. The mean value (S) is parallel to the Sy axis. 

The minimum- variance determination prescribed by P2 is 
also an optimal strategy for polarization- squeezing detec- 
tion iEHi . In this case, it suffices to consider a general Stokes 
parameter rotated by in the dark plane (orthogonal to the 
direction of (S)), namely, .§61 = cos6 Sx + sin Sz, so that 
(Se) = 0. In addition, since for bright fields the fluctu- 
ations are small compared with the mean values, one has 
{ASe)^ - lN[{AXH,e)^ + {AXv,e)% where Xg are the 
rotated quadratures for each polarization mode. The searched 
point is obtained by optimizing over 0, finding 2° ± 0.3° 1I22I] . 

From the data we get Pi = 1 and P2 ^ 1 (within the exper- 
imental precision). This is simply due to the large excitation 
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of the Stokes vector, which dominates in the definition (|5]). 
In the next experiment, we produce a state with a very small 
excitation, so that the second-order degree is governed by the 
Stokes fluctuations. 

We generate a quadrature squeezed vacuum in a well- 
defined spatio-temporal polarization mode using an OPO op- 
erating below threshold and pumped with a 532 nm light beam 
(see Fig 2) 1I23I1 . The parametric down-conversion interaction 
is based on a type I phase-matched periodically poled KTP 
crystal, which produces a squeezed vacuum in the H mode 
while leaving the V mode in the vacuum. With this informa- 
tion at hand, we write the resulting state as qh ^ \0)v v{0\ 
where qh is the density operator of the state produced in the 
OPO. 

In contrast to the previous experiment, now we characterize 
the polarization by using two-mode homodyne detection. As 
this provides complete knowledge about the measured state, 
the Stokes fluctuations will be contained in the homodyne 
data. Since the H and V modes are known to be uncorre- 
cted, a complete reconstruction can be obtained just by mea- 
suring the two orthogonal modes independently. To this end, 
we direct our two-mode state to a standard homodyne detec- 
tor where the polarization of the local oscillator (LO) could be 
swapped between H and V polarizations. The measurements 
results are demodulated at a sideband frequency of 5 MHz 
with a bandwidth of 100 kHz. The total detection efficiency is 
about 87 %. 

Using the time-resolved data of the orthogonal modes as 
well as the a priori state information, we fully reconstruct 
the density matrix in a 16-dimensional Fock space using a 
maximum likelihood algorithm. From this density matrix we 
calculate the moments of the Stokes parameters and plot 
the result as a color map on the Poincare sphere, as shown in 
Fig 2. Now N ^ l.b and the degree of polarization of the 
quadrature squeezed vacuum state is not only governed by the 
first moment (as the bright squeezed state). This is nicely il- 
lustrated in the new definition of the degree of polarization: 
we calculate Pi = 0.998 ± 0.001 and P2 = 0.79 ± 0.01. We 
also note that P2 = 0.43 ± 0.01 and it does not capture the 
variations on the Poincare sphere. 

Discussion and concluding remarks. — The definition (O 
has proved to be a satisfactory solution to deal with second- 
order polarization properties. Of course, a complete character- 
ization must involve a whole hierarchy of polarization degrees 
P/c containing all the orders, as it happens with field correla- 
tions in coherence theory. Although the second order consid- 
ered here surely accounts for most of the interesting, and in 
many cases dominant effects, some subtleties may arise when 
taking into account higher orders. A full analysis of these 
questions exceeds the scope of this work and will be presented 
elsewhere. 

Our approach is based on the underlying SU(2) symmetry 
of light polarization. This makes possible a direct transla- 
tion of our results to other fields where the same symmetry 



plays an important role, such as cold atoms i24ll . The ap- 
proach is also well suited for other unitary symmetries, such 
as SU(2)*^^ or SU(3). The former is connected with the polar- 
ization of spatial-multimode fields [5], while the latter has re- 
cently attracted a lot of attention in relation with near-field op- 
tics 1I25I1 . In these cases, the optimization process can be more 
involved, but the spirit of our approach remains the same. 
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